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ABSTRACT 

Context. Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered 
by young and energetic pulsars. HESS 11857+026 is a spatially extended y- ray source detected by H.E.S.S. and classified as a possible 
pulsar wind nebula candidate powered by PSR 11856+0245. 

Aims. We search for y-ray pulsations from PSR 1 1 856+0245 and explore the characteristics of its associated pulsar wind nebula. 
Methods. Using a rotational ephemeris obtained from the Lovell telescope at lodrell Bank Observatory at 1.5 GHz, we phase-fold 
36 months of y-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete y-ray spectral and 
morphological analysis. 

Results. No pulsation was detected from PSR J1856+0245. However, significant emission is detected at a position coincident with 
the TeV source HESS 11857+026. The y-ray spectrum is well described by a simple power law with a spectral index of F = 1.53 ± 
0-11 stat ± 0.55 syst and an energy flux of G(0.1-100 GeV) = (2.71 ± 0.52 stat + 1 .5 l sy st) x 10 -11 ergs cm -2 s _1 . This implies a y-ray 

efficiency of ~ 5 %, assuming a distance of 9 kpc, the y-ray luminosity of Lj/ WJV (0. 1-100 GeV) = (2.5 ± 0.5 sm , ± 1 x 10 35 ( 9 ^:) 
ergs s _1 and E = 4.6 x 10 36 erg s -1 , in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new 
constraints on its pulsar wind nebula nature. 

Key words, pulsars : general, pulsars : individual object : PSR 11856+0245, ISM : individual object : HESS 11857+026, Gamma 
rays : general 


1. Introduction 


Pulsar wind nebulae (PWNe) are bubbles of shocked relativistic 
particles produced by the interaction of the pulsar’s wind with 
the surrounding medium (Gaensler el al. 2006). Since 2003, the 
continuous observations of the Galactic Plane by Cerenkov tele- 
scopes have yielded the detection of more than 60 Galactic TeV 
sources. Among them, PWNe are the dominant class with 29 
firm identifications. In the GeV energy range, 5 PWNe have been 
firmly identified by the Fermi- LAT. They are all powered by en- 


ergetic pulsars and their y-efficiencies are ~ 1%, consistent with 
TeV observations (Ackermann et al. 2011). 

The presence of a pulsar close to the source position is an 
important clue to confirm the identification of a PWN, which 
often requires information from the radio/X-ray wavelengths. 
Radio/X-ray PWNe are often associated with TeV extended 
sources offset from their pulsars, which can be explained by 
an inhomogeneous environment (Hinton et al. 2010). In such 
sources, TeV radiation can be explained by Inverse Compton 
(IC) scattering of accelerated leptons on ambient photon fields 
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(CMB, IR, ...) or by n <] decay from the interaction of accelerated 
hadrons with nuclei of the interstellar medium. 

HESS J1857+026 is a very high energy (VHE) y-ray 
source detected by H.E.S.S. during the Galactic Plane Survey 
(Aharonian et al. 2008). The extended (~ 0.11°) TeV source 
was identified as a PWN candidate after the discovery of PSR 
J1856+0245 (offset ~ 0.12°) in the Arecibo PALFA survey 
(Hessels et al. 2008) with a dispersion measure of 222 cm -3 
pc. Recently, MAGIC reported a measured extension in the 0.2- 
1 TeV energy range significantly larger (0.22°) than the ex- 
tension reported by H.E.S.S. in the 0.6-80 TeV energy range 
(Klepser et al. 2011). PSR J1856+0245 is an energetic pulsar 
(E = 4.6 x 10 36 erg j _1 ) located in a crowded region, 1.3° from 
the bright SNR W44 (Abdo et al. 2010) and 0.6° from the fainter 
SNR HESS J1 858+020 on which only an upper limit could be 
set using Fermi data (Torres et al. 2011). Significant emission 
coincident with HESS J 1857+026 was observed above 100 GeV 
using Fermi-LAI observations (Neronov et al. 2010). 

Here, we report in detail GeV observations of the HESS 
J1857+026/PSR J1856+0245 system using Fermi-LAI observa- 
tions and discuss their implications for the nature of the source. 

2. LAT description and data selection 

The LAT is a y-ray telescope that detects photons by conversion 
into electron-positron pairs and operates in the energy range be- 
tween 20 MeV and 300 GeV. Details of the instrument and data 
processing are given in Atwood et al. (2009). The on-orbit cali- 
bration is described in Abdo et al. (2009a). 

The following analysis was performed using 36 months of 
data collected from August 4, 2008 to August 31, 201 1 within a 
10 x 10° square around the position of HESS J 1857+026 aligned 
with Galactic coordinates. We excluded y-rays coming from a 
zenith angle larger than 100° because of possible contamination 
from secondary y-rays from the Earth’s atmosphere (Abdo et 
al. 2009b). We used the P7 V6 Instrument Response Functions 
(IRFs), and selected the ‘Source’ events which correspond to the 
best compromise between the number of selected photons and 
the charged particle residual background for the study of point- 
like or slightly extended sources. 

3. Data analysis 

3. 1. Timing analysis of PSR J 1856+0245 

With its large spin-down power, PSR J 1 856+0245 is one of the 
more energetic radio pulsars known. Its spin period of 80.9 ms 
and characteristic age of 20.6 kyr are similar to those of the 
Vela pulsar. The dispersion measure and NE2001 electron den- 
sity model of the Galaxy assign PSR J 1856+0245 a distance of 
~ 9 kpc (Cordes et al. 2002). 

This pulsar is not monitored as part of the LAT pulsar timing 
campaign (Smith et al. 2008), as it was discovered subsequently, 
but has nevertheless been regularly observed with the Lovell 
telescope at Jodrell Bank Observatory (Hobbs et al. 2004). The 
ephemeris of PSR J 1 856+0245 used in the analysis was obtained 
using 82 observations at 1.5 GHz made with the Lovell telescope 
between May 4, 2008 and September 4, 2011. The arrival times 
of events were corrected to the Solar System Barycenter using 
the JPL DE405 Solar System ephemeris. The TEMP02 timing 
package (Hobbs et al. 2006) was then used to build the timing 
solution. We fit the radio times of arrival (TOAs) to the pulsar 
rotation frequency and first four derivatives (in order to remove 
timing noise). By including a fourth derivative the RMS of the 


Parameter 

Value 

v (Hz) 

12.3597551142(1) 

v (1) (Hz s _1 ) 

-9.48698(1) x 10“ 12 

v (2) (Hz s -2 ) 

1.6585(9) x 10“ 22 

v (3) (Hz s -3 ) 

2.47(7) x 10- 30 

v (4) (Hz s“ 4 ) 

1.32(7) x 10- 37 

DM (cm -3 pc) 

622 

Period epoch (MJD) 

55128 

Start time (MJD) 

54615 

End time (MJD) 

54570 

Number of TOAs 

82 

TOA rms (ms) 

1.08 


Table 1 . Parameters of the fit of the TOAs. v correspond to the rotational 
frequency of the pulsar and v (,) its time derivative of order i. In parenthe- 
ses are the 1 cr uncertainty on the least-significant digits quoted, from 
TEMP02. 


timing residuals decreases in about 30%. It goes from 1.551 ms 
to 1.081 ms, which is equivalent to a decrease from 19 to 13 mil- 
liperiods. The reduced^ 2 of the fit also decreases significantly, 
changing from 22 to 5. The fit results are summarized in Table 

1 . 

This timing solution will be made available through the 
Fermi Science Support Center 1 (FSSC). 

For the LAT analysis, photons with energies above 100 
MeV and within a radius of 1.0° of the radio pulsar posi- 
tion <+(12000) = 18 A 56'"50.937 I , <5(J2000) = +02°45' 47.046" 
were selected using an energy-dependent cone of radius 6 < 
max(5.12° x (.E/lOOMeV) - , 0.2°) and phase-folded using the 
radio ephemeris previously described. This choice takes into ac- 
count the instrument performance and improves the signal to 
noise ratio over a broad energy range. No significant pulsation 
was detected for all tested energy bands (100 MeV - 300 GeV, 
100 MeV - 300 MeV, 300 MeV - 1 GeV, > 1 GeV). Following 
the procedure used by Romani et al. (2011), we fitted a point 
source at the position of PSR J1 856+0245 in the 0.1-1 GeV 
energy range assuming a power-law of index 1.62 and a cut-off 
energy at 2.8 GeV to derive a 99% Bayesian upper limit on the 
flux of 3.27 x 10 -8 ph cm -2 s _1 , well below typical y- ray fluxes 
reported for pulsars detected by Fermi-LAI (Abdo et al. 2010c). 


3.2. Spatial and spectral analysis 

Two different tools were used to perform the spatial and spec- 
tral analysis: gtlike and pointlike, gtlike is a binned 
maximum-likelihood method (Mattox et al. 1996) implemented 
in the Science Tools distributed by the FSSC. pointlike is 
an alternate binned likelihood technique, optimized for char- 
acterizing the extension of a source (unlike gtlike), that was 
extensively tested against gtlike (Kerr 2011). These tools fit 
a source model to the data along with models for the instru- 
mental, extragalactic and Galactic components of the back- 
ground. In the following analysis, the Galactic diffuse emission 
is modeled by the standard LAT diffuse emission ring-hybrid 
model gaL2yearp7v6_v0.fits. The residual cosmic-ray back- 
ground and extragalactic radiation are described by a single 
isotropic component with a spectral shape described by the 
file iso_p7v6source.txt. The models have been released and de- 
scribed by the Fermi-LAI Collaboration through the FSSC 2 . 

1 FSSC:http://fermi. gsfc.nasa.gov/ssc/data/access/lat/ephems/ 

2 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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The source significance is measured by a test statistic (TS) 
defined as TS= 2(log(Li) - log(L 0 )), where L\ corresponds to 
the likelihood obtained by fitting a model of the source of inter- 
est and the background model and Lq corresponds to the likeli- 
hood obtained by fitting the background model only. In the fol- 
lowing, the correspondence between the significance and the TS 
value is evaluated from the^ 2 distribution with 4 degrees of free- 
dom (position and spectral parameters). 

The 41 sources within 15° of HESS J 1857+026 in the 
Second Fermi - LAT catalog (Abdo et al. 2011) were taken into 
account. We refitted all spectral parameters of the 16 sources 
within 5° around HESS J 1857+026. The region includes the 
bright SNR W44, known to interact with its environment. 
Extended and only 1.3° from HESS 11857+026, W44 could in- 
fluence our fit. We refitted it assuming an elliptical ring and 
obtained results consistent with those of Abdo et al. (2010a). 
The centroid is located at a = \S h 56 m , 5 = +01°22'. The fit- 
ted semi major and semi minor axes are respectively maj/2 = 
(0.33 ± 0.10 s(ar )°, min/2 = (0.20 ± 0.02 stor ) 0 . The angle of the 
semi major axis from celestial North, considered positive toward 
increasing right ascension, is (327 ± 22 stat )° . 

3.2.1. Shape and position of HESS J1 857+026 counterpart 

Source shape analysis requires the best possible angular resolu- 
tion. Since the source has a hard spectrum (see Section 3.2.2) we 
made a compromise between statistics and resolution by select- 
ing photons above 10 GeV. This drastically reduces the contribu- 
tion of the Galactic diffuse background and improve the single- 
photon angular resolution. Fig. 1 (Top) presents a LAT TS Map 
in the energy range of 10 GeV to 300 GeV. To each pixel is 
associated a TS value calculated assuming a point source in its 
center and fitting only the flux of the source assuming a power 
law spectrum with a spectral index of 2. A source coincident 
with HESS J1857+026 is clearly visible. We determined the ex- 
tension of the source using point like with three different mod- 
els : a point source, a uniform disk and a Gaussian. No signifi- 
cant extension was obtained above 10 GeV. The GeV emission 
was fit to position n(72000) = 18''57 m , (5(72000) = +02°45' 
with an average statistical error of 0.05°, consistent with the 
position determined by H.E.S.S. ff(72000) = 18 A 56"'50.80" , 
5(72000) = +02°45'50.2". 

3.2.2. Spectral analysis 

Fig. 1 (Middle) shows a TS map of the region in the energy range 
0.1-1. 3 GeV. There is excess emission near HESS J1857+026 
located at a(72000) = 18 ft 54 m , 5(72000) = +02°59'. This ex- 
cess is inconsistent with that of HESS J 1 857+026 and was added 
to the background model. This additional background source 
was fitted assuming a pure power-law with an integrated flux 
of F(0.1-100 GeV)=(2.31 ± 0.37 stat ) x 10“ 7 ph cm” 2 s“\ a spec- 
tral index of T = 3. 1 8 ± 0.56 stat , which gives a significance above 
300 MeV of ~ 3.6cr (TS=20). 

Spectral analysis was performed using gtlike, selecting 
only 0.3-300 GeV to avoid the low energy range that is dom- 
inated by the diffuse Galactic background and subject to high 
systematics. In this energy range, HESS J 1857+026 is well de- 
scribed by a pure power-law with an integrated flux extrapo- 
lated down to 100 MeV of F(0. 1-100 GeV)=(5.78 ± 0.85 stat ± 
3.11 syst) x 10 -9 ph cm -2 s -1 , a spectral index of T = 1.53 ± 
0.1 l stat ±0.55 syst and an energy flux of G(0. 1-100 GeV)=(2.71 ± 
0.52 stat ± 1 .5 l S yst) x 10 -11 ergs cm -2 s _1 , which gives a signifi- 


cance above 300 MeV of ~ 5.4<x (TS=39). The residual TS map 
after fit is presented in Fig.l (Bottom) and shows no significant 
excess. 

Fermi - LAT spectral points for HESS J1857+026 were ob- 
tained by splitting the 0.3-100 GeV range into 4 logarithmically- 
spaced energy bins plus a bin between 100 and 300 GeV which 
contains 22 photons corresponding to a TS of 14, as presented in 
Fig. 2. A 99 % C.L. upper limit is computed when TS<10 using 
the approach used by Abdo et al. (201 1). The errors on the spec- 
tral points represent the statistical and systematic uncertainties 
added in quadrature. 

Three main systematic uncertainties can affect the LAT flux 
estimate for a point source: uncertainties in the Galactic diffuse 
background, uncertainties on the effective area and uncertain- 
ties on the shape of the source. The dominant uncertainty at low 
energy comes from the Galactic diffuse emission, estimated by 
changing the normalization of the Galactic diffuse model artifi- 
cially by ±6% as done in (Abdo 2010 b). Since it is computed for 
P6 IRFs, the 6% factor overestimates the bias. The second sys- 
tematic is estimated by using modified IRFs. The fact that we do 
not know the true y-ray morphology introduces a last source of 
error. We derived an estimate of the uncertainty on the shape of 
the source by using the best Gaussian model obtained by HESS. 
We combine this various errors in quadrature to obtain our best 
estimate of the total systematic error at each energy and propa- 
gate through to the fit model parameters. 

Assuming a distance of 9 kpc and an isotropic emission, the 
y-ray flux corresponds to L y pwN (S). 1-100 GeV) = (2.5 ± 0.5 Jta/ ± 

1.5 jyjf ) x lO 35 ^)- ergs s 1 . Using the pulsar’s E, this yields 
a y-ray efficiency of ~ 5 %, which is one of the highest PWN 
efficiencies observed at GeV energies (Ackermann et al. 2011). 
It is still in the range of expected values for PWNe seen by Fermi 
and is close to the estimate of 3% using HESS data (Mattana et 
al. 2009). 


4. Supporting X-Ray measurement 

To obtain a precise flux for any potential X-ray PWN associated 
with PSR 11856+0245, we analyzed a 39-ks Chandra ACIS-I 
observation from February 28, 2011 (Obs. ID 12557). These data 
were recorded in the VFAINT and Timed Exposure (TE) modes 
and were analyzed using CIAO 3 * version 4.3.1 with CALDB 
4.4.3. PSR J1856+0245 is clearly detected as a point source, 
but there was no immediate evidence for extended emission sur- 
rounding this position. In the 2-10 keV energy range, assuming 
a power law spectrum to fit the XMM data, the unabsorbed flux 
of the pulsar is 8.3+ 2 g x 10 -14 erg s -1 . An in-depth analysis of 
the X-ray properties of PSR J 1 856+0245 as well as the XMM 
data will be presented in Bogdanov et al. (in preparation). 

Given that the size of the potential X-ray PWN is not known, 
we investigated an extraction region to see whether they produce 
a statistically significant excess of counts compared with the 
background. This extraction region is in the form of annuli ex- 
tending from 2" - 15" from the position of the pulsar. The inner 
radius was chosen to avoid contamination from the pulsar. The 
outer radius was chosen based on the X-Ray PWNe observed 
for pulsars with comparable E by scaling their angular size with 
their distance (Kargaltsev & Pavlov 2008). The background re- 
gions were chosen from several other source-free regions in the 
vicinity of the pulsar. For the 2" - 15" extraction region we 

3 Chandra Interactive Analysis of Observations (Fruscione et al. 

2006). 
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find an upper limit on the unabsorbed flux of 5 x 10“ 14 erg s -1 
cm -2 (l-10keV, 3 <x confidence), corresponding to a luminosity 
of 5 x 10 32 erg s -1 . The counts in this region show a 2cr ex- 
cess from zero counts. Assuming now the TeV position and a 
6' extraction region consistent with the HESS morphology, we 
find an upper limit on the unabsorbed flux of 2 x 10 -12 erg s -1 
cm -2 (l-10keV, 3 cr confidence), corresponding to a luminosity 
of 2 x 10 34 erg s _1 . Given the marginal significance of the count 
excesses derived, we cannot convincingly claim the detection of 
a weak X-ray PWN. These luminosity limits are derived from 
the 3cr upper bound on the net count rate and assume a typical 
power-law spectrum of index 1 .5 for the PWN, a distance of 9 
kpc, and a column density V H = 4 x 10 22 cm -2 based on the 
spectroscopic analysis of the XMM data. 

5. Discussion 

To investigate the global properties of the PWN, we apply a 
one-zone time dependent SED model which reproduces the 
multi-wavelength measurements from MAGIC (Klepser et 
al. 2011), H.E.S.S. (Aharonian et al. 2008), as well as the LAT 
and the 2" - 15" X-ray upper limit described above. This 
model is described in Grondin et al. (2011) and Abdo et al. 
( 2010 ), with a more detailed description (applied to a multi- 
zone scenario rather than one-zone) found in Section 3.10 of 
Van Etten & Romani (2011). 

The model computes SEDs from evolving electron pop- 
ulations over the lifetime of the pulsar in a series of time 
steps with the energy content of the injected particle popula- 
tion varying with time according to the pulsar spin down. 
At each time step, new particles are injected and syn- 
chrotron, inverse Compton, and adiabatic cooling is com- 
puted. Synchrotron and inverse Compton fluxes are com- 
puted from the final electron spectrum, magnetic held, 
and photon fields. Starting assumptions include: the source 
of all particles and magnetic field in the nebula is PSR 
J1856+0254, ambient photon fields are static, and the neb- 
ula is 9 kpc distant. 

The large size of the nebula (~ 20 pc for a distance of 
9 kpc) implies a middle-aged PWN, whose parent SNR has 
likely evolved from the free-expansion phase to the Sedov 
phase. The Sedov phase is expected to occur on a timescale of 
tsedov ~ 7 kyr for an explosion of 10 51 erg, an ejecta mass of 
10MO, and an ambient medium density of 1 cm -3 (Reynolds 
& Chevalier 1984). Eventually, the inward moving SNR re- 
verse shock collides with the expanding PWN, which can 
happen as late as late as 5 times the transition to the Sedov 
phase (van der Swaluw et al. 2004). The interaction of the 
PWN and the SNR reverse shock compresses the PWN, re- 
sulting in an increased magnetic held. Even in the spherically 
symmetric case the evolution is complex, with the PWN un- 
dergoing a series of oscillations due to the reverse shock in- 
teraction. The significant offset of the pulsar from the y-rav 
centroid implies either a PWN expansion into an inhomoge- 
nous medium, or an asymmetric reverse shock interaction. 
Either way, models such as Gelfand et al. (2009) which com- 
pute the PWN compression assuming spherical symmetry 
cannot be applied. Given the lack of multiwavelength data 
for HESS J 1857 +026 and the complexity in SNR-PWN in- 
teraction, we see little reason to adopt a complicated spa- 
tial model and therefore assume a simple spatial evolution. 
During the free-expansion phase of the PWN (assumed to 
be ~ 10 4 years) we adopt an expansion of R oc t, following 
which the radius evolves as R oc t° \ appropriate for a PWN 


expanding in pressure equilibrium with a Sedov phase SNR 
(Reynolds & Chevalier 1984). We also adopt a simple mag- 
netic field model, and assume that over the pulsar lifetime 
the magnetic field evolves as B oc f -1 - 5 , following ~ 500 years 
of constancy. This magnetic field evolution is similar to the 
r 1 3 - t~ 2 behavior derived by Reynolds & Chevalier (1984) 
and the f -1 7 evolution of Gelfand et al. (2009) during the ini- 
tial expansion phase. The power-law evolution of magnetic 
field is similar to that adopted in Zhang et al. (2008), Lemiere 
et al. (2009), Mayer et al. (2012) for similarly aged PWNe. 

Model fitting is achieved by minimizing the y 2 between 
model and data using the downhill simplex method described 
in Press et al. (1992). For each ensemble of N variable pa- 
rameters we evolve the system over the pulsar lifetime and 
calculate y 2 between model curves and flux data points. The 
simplex routine subsequently varies the parameters of inter- 
est to minimize the fit statistic. We estimate parameter er- 
rors by computing x 2 for a sampling of points near the best 
fit values and using these points to fit the A- dimensional 
ellipsoid describing the surface of Ay 2 = 2.71. Assuming 
that the adopted PWN model is correct, and that errors are 
Gaussian, the projected size of this Ay 2 = 2.71 ellipsoid onto 
each parameter axis defines the 90% multi-parameter (pro- 
jected) error. 

Using our ephemeris, we tried to derive a value for the brak- 
ing index using eq. 1 where v, v® and n represents respectively 
the rotational frequency of the pulsar, its time derivative of order 
i and the braking index. 


We obtained n>20. Large braking indices between glitches 
are common among Vela-like pulsars and are likely to be asso- 
ciated with glitch recoveries. These large values should not be 
interpreted as the long-term braking index due to secular spin 
evolution but instead these correspond to transient states caused 
by large glitch activity as discussed in section 3.2.2 of Hobbs 
et al. (2010). Dipole braking indices have been measured only 
for a few pulsars with the highest spindown rates (see Table 1 of 
Espinoza et al. (201 1)). 

With its ~21 kyr of characteristic age, PSR J1856+0245 is a 
Vela-like pulsar possibly affected by quite a large glitch activity, 
which we have not yet directly seen. Thus, we fix here the pulsar 
braking index to the canonical value of 3. 

We assume the existence of three primary photon fields 
(CMBR, far IR (dust), and starlight) and use the interstel- 
lar radiation held from GALPROP (Porter et al. 2005) to 
estimate the photon fields at the Galactic radius of PSR 
J1856+0245. A distance of 9 kpc in the direction of the pul- 
sar corresponds to a Galactic radius of 5.4 kpc. At this ra- 
dius, the peak of the SED of dust IR photons corresponds to 
a black body temperature of T ~ 32 K with a density of ~ 1.1 
eV cm -3 , while the SED of stellar photons peaks at T ~ 2500 
K with a density of ~ 1.2 eV cm -3 . 

A simple exponentially cutoff power-law injection of elec- 
trons, evolved properly over the pulsar lifetime, often pro- 
vides an adequate match to PWN SEDs. For this injection 
spectrum we fit four variables: final magnetic field B f , elec- 
tron high energy cutoff £ cut , electron power-law index p, and 
initial pulsar spin period Pq ms. The best fit parameters, with 
errors, are given in Table 2. The initial spin period gives an 

3 http://galprop.stanford.edu/resources.php?option=data 
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age of 20 ± 1 kyr. This model poorly matches the low energy 
MAGIC points, as shown in Figure 2 (Top). 

Another option to fit the multi-wavelength data is to 
adopt the relativistic Maxwellian plus power-law tail elec- 
tron spectrum proposed by Spitkovsky (2008). We imple- 
ment this spectrum as described in Grondin et al. (2011). The 
best fit, presented in Fig. 2 (Middle) and Table 2, is obtained 
with kT = 0.64 TeV corresponding to an upstream Lorentz 
factor of 2.5 x 10 6 . A very high cutoff of 450 TeV is required, 
though the power-law index of p = 2.49 is consistent with the 
value of ~ 2.5 proposed by Spitkovsky (2008). The initial spin 
period gives an age of 13 ± 1 kyr. The relativistic Maxwellian 
plus power law model better matches the multi-wavelength 
data, and also directly probes the upstream pulsar wind via 
fitting of the upstream Lorentz factor of the wind. 

A hadronic scenario is also possible, with y-rays aris- 
ing from proton-proton interactions. For this model, corre- 
sponding to Fig. 2 (bottom) and Table 2, we fix the ambient 
gas density at 50 cm -3 and age at 20 kyr. 

6. Conclusions 

Using 3 years of Fermi-LAI data, a y-ray source has been de- 
tected at high significance at a position coincident with the PWN 
candidate HESS J1857+026. This VHE spectrum observed by 
MAGIC and H.E.S.S., combined with the limits imposed by the 
steep LAT data, is difficult to match with a simple power-law in- 
jection of electrons (or protons), and we find a significantly bet- 
ter fit with a relativistic Maxwellian plus power-law spectrum. 
The low magnetic field of the leptonic fits, due to the stringent X- 
ray upper limit, implies that if PWN leptons are indeed responsi- 
ble for the y-ray flux, they must be dominated by relic electrons 
which have escaped the PWN core into weakly magnetized sur- 
roundings. The hadronic scenario relaxes this constraint, though 
the energy requirements are quite high even for a dense ambi- 
ent medium, and a very hard power-law index is required. At 
present the true nature of HESS J 1857+026 remains a mystery, 
though the new LAT data and X-ray upper limit hint that this 
source may be another relic PWN, increasing the population of 
such high energy y-ray systems. 
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Fig. 1. TS maps computed by pointlike. The green crosses represent 
the sources of the 2FGL catalog included in the model, whereas the 
blue X represents the source we added in the model. The green contours 
represent the H.E.S.S. data (Aharonian et al., 2008). The magenta cir- 
cle represents the position of PSR J1856+0245. Top: TS map obtained 
between 10 and 300 GeV. The position of the Fermi excess is consis- 
tent with that of H.E.S.S. Note that HESS J1857+026 is not included 
in the model. Middle: TS map obtained between 0.1 and 1.3 GeV. 
This figure shows the residual excess taken into account in our model. 
Bottom: Residual TS map obtained between 10 and 300 GeV when all 
the sources are taken into account. 
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Fig. 2. Spectral energy distribution of HESS J1857+026 with a simple 
exponentially cutoff power-law electron spectrum (Top), a relativistic 
Maxwellian plus power-law electron spectrum (Middle), and an expo- 
nentially cutoff power-law proton spectrum (Bottom). The X-ray flux 
upper limit obtained using Chandra(green), LAT spectral points (red), 
MAGIC points (violet) (Klepser et al. 2011), and H.E.S.S. points (blue) 
(Aharonian et al. 2008) are shown. The black line denotes the total 
synchrotron, inverse Compton and pion decay emission from the neb- 
ula.Thin curves indicate the Compton components from scattering on 
the CMB (long-dashed), IR (medium-dashed), and stellar (dotted) pho- 
tons. 
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Table 2. Model Parameters 


Model 

n a 

edov (k yr) 

/3(B oc /) 

B f (n G) 

P (E~n 

kT (TeV) / n(cm _3 ) c 

£cut (TeV) 

Po(ms)/E 0 (xl0 50 erg) d 

^/d.o.f. 

Power-Law 

y 

10' 

-1.5' 

3.7 ± 0.5 

2.08 ± 0.03 

- 

99 + 37 

9.6 ±6.1 

24.9/21 

Rel. Max. 

y 

10' 

-1.5' 

2.6 ± 0.7 

2.49 ± 0.09 

0.64 ± 0.09 

450 ± 290 

36 ±4 

12.4/20 

Hadron 

y 

10' 

-1.5' 

20+8° 

z,u -20 

1.83 ±0.04 

50“ 

75 + 25 

0.64 ± 0.06 rf 

24.6/21 


“ Pulsar braking index 
b Sedov phase onset, after which r oc t 03 
1 Ambient medium density 
d Total energy injected 
' Held fixed 
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